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"^ ABSTRACT 

(N 

Context. Planets have been detected around a variety of stars, including low-mass objects, such as brown dwarfs. However, such 
extreme cases are challenging for planet formation models. Recent sub-millimeter observations of disks around brown dwarf mea- 
sured low spectral indices of the continuum emission that suggest that dust grains grow to mm-sizes even in these very low mass 
environments. 

Aims. To understand the first steps of planet formation in scaled-down versions of T-Tauri disks, we investigate the physical conditions 
that can theoretically explain the growth from interstellar dust to millimeter-sized grains in disks around brown dwarf. 
Methods. We modeled the evolution of dust particles under conditions of low-mass disks around brown dwarfs. We used coagulation, 
Q fragmentation and disk-structure models to simulate the evolution of dust, with zero and non-zero radial drift. For the non-zero radial 

?-H drift, we considered strong inhomogeneities in the gas surface density profile that mimic long-lived pressure bumps in the disk. We 

jyFj studied different scenarios that could lead to an agreement between theoretical models and the spectral slope found by millimeter 

Ci observations. 

' ' Results. We find that fragmentation is less likely and rapid inward drift is more significant for particles in brown dwarf disks than 

in T-Tauri disks. We present different scenarios that can nevertheless explain millimeter-sized grains. As an example, a model that 
'"^ combines the following parameters can fit the millimeter fluxes measured for brown dwarf disks: strong pressure inhomogeneities of 

^ ~ 40% of amplitude, a small radial extent ~ 15 AU, a moderate turbulence strength ttiuib = 10"^, and average fragmentation velocities 

^^^ for ices I'f = 10 mi-"'. 

cn ^ 

\^ Key words, accretion, accretion disk - circumstellar matter -(stars:) brown dwarf-planet formation. 

f^ 1 . Introduction 2009), a flat tendency for the disk geometry inferred from Spitzer 

C^ observations and spectral energy distribution (SED) modeling 

Since the first confirmed discovery of brown dwarf (BD)re/£/e 7 (e.g. Apai et al. 2005; Allers et al. 2006; Guieu et al. 2007 

(Rebolo et al. 1995) and Gliese 229B (Nakajima et al. 1995) Scholz et al. 2007; Morrow et al. 2008; Pascucci et al. 2009 

in 1995, several hundred BDs have been identified and many Szucs et al. 2010), and a longer lifetime (Carpenter et al. 2006 

efforts have focused on understanding these objects which are rj^z et al. 2012; Harvey et al. 2012b). 
in considered to be an intermediate step between planets and stars. Xo study dust growth in protoplanetary disks, different mech- 

C^ Observations of BD show near-infrared excess emission (e.g. anisms should be taken into account, as for instance turbulent 

Muench et al. 2001; Liu et al. 2003) that reveals material around mixing, settling, aerodynamical drag with the gas, collision rates 

young BDs. Moreover, typical fluxes measured with millimeter and fragmentation. Essentially, when the dust particles are smaU, 

observations (e.g Klein et al. 2003; Scholz et al. 2006; Joergens they are well coupled to the gas, move along with it, and grow 

et al. 2012) are in most of the cases lower than few mJy at 1 mm, as a consequence of surface forces. However, when the particles 

implying that these circumsteUar disks have masses of few Mjup increase in size, they start to decouple from the gas and the rel- 

or even lower. Determining whether these low-mass disks can ative velocities between particles increase, leading to fragmen- 

be the scene of the formation of planetesimals or even planets is tation collisions (Weidenschilling 1977; Brauer et al. 2008). In 

still a subject of debate. addition, before any meter-sized object can be formed at Earth- 

Observationally, some aspects of disks around BD are dif- Sun distances, dust drifts toward the central star because of the 

ferent from those around T-Tauri and Herbig Ae/Be stars. They sub-Keplerian velocity of the gas. Millimeter (mm) grains expe- 

have a lower accretion rate (~ 10"'^ Mq yr"', Herczeg et al. rience the same rapid inward migration in the outer regions of 
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the disks, even though mm-size particles have been observed in 
those regions of protoplanetary disks (e.g. Wilner et al. 2000; 
Natta et al. 2004; Rodmann et al. 2006; Ricci et al. 2010a,b, 
2011; Ubach et al. 2012). Local pressure maxima in disks have 
been proposed as a solution of this rapid inward drift (Klahr & 
Henning 1997; Johansen et al. 2009; Pinilla et al. 2012a). 

Birnstiel et al. (2010b) showed that under typical T-Tauri pa- 
rameters, neglecting radial drift and considering different param- 
eters for the disk such as turbulence, particles could reach mil- 
limeter sizes, meaning that the spectral index Q'i-3mm can have 
low values. Birnstiel et al. (2010b) models predict that for a disk 
mass lower than 5 Mjup, such as the mass of a BD disk, a'i_3n,jn 
would be close to 3. However, recent millimeter observations 
confirmed low values ai-3mm for two BD disks (Ricci et al. 
2012; Ricci et al. 2013), with ai_3mni ~ 2.3. In addition, as we 
show below, the radial drift and fragmentation barriers are differ- 
ent for particles in disks around BD than for those around more 
massive and luminous stars. Radial drift indeed has a stronger 
influence on particles in BD disks, and as a consequence, any 
mechanism in a disk that may allow a reduced inward migra- 
tion of grains has to be more extreme in BD disks to lead to an 
effective trapping of particles. Explaining how the first pebbles 
are formed from interstellar dust in BD disks is therefore a very 
intriguing topic. 

The purpose of this paper is to investigate whether the dust 
growth models that were successful for T-Tauri disks are, when 
applied to BD disks, consistent with a short set of observations at 
millimeter wavelengths (two measurements from recent ALMA 
and CARMA observations). As mentioned above, infrared stud- 
ies have on the other hand focused on larger surveys, showing 
interesting trends such as a flatter BD disk geometry. However, 
we do not aim to reproduce these trends, because infrared obser- 
vations yield no strong constraints on the dust evolution models 
(other than those studied in e.g. Szucs et al. (2010) or Mulders 
&Dominik(2012)). 

For this work, dust coagulation/fragmentation models were 
considered in two main cases: first, the extreme case of setting 
the radial drift to zero, and second, taking radial drift into ac- 
count with strong inhomogeneities in the gas surface density 
that mimic long-lived pressure bumps. In both cases, different 
disk parameters were considered, to analyze which scenarios 
are the best incubators of the first pebbles found in BD disks 
with millimeter observations. In Sect. 2, we describe the drift 
and fragmentation barriers for the specific case of BD disks and 
the physical parameters of the dust coagulation/fragmentation 
model. Numerical results, observational perspectives, and com- 
parisons with current mm-observations are presented in Sect 3. 
A summary of the results with the corresponding discussion of 
this work is recapitulated in Sect. 4. Finally, the main conclusion 
of this paper is given in Sect. 5. 



2. Dust evolution model, drift, and fragmentation 
barriers in BD disks 



relative velocities substantially increase through turbulent mo- 
tion and radial drift; the collision energies are therefore high 
enough to cause fragmentation, as was experimentally shown by 
Blum & Wurm (2008). From the theoretical point of view, stellar 
and disk properties strongly influence these first steps of planet 
formation, meaning that effects such as the usual inward migra- 
tion of dust grains may vary for BD disks. In this section, we 
first explain the main characteristics of the numerical model for 
the dust evolution, followed by an explanation of the drift and 
fragmentation barriers in BD disks, and finally we describe the 
set-up for the numerical simulations that are considered in this 
work. 



2.1. Dust evolution model 

For the dust evolution, we used the coagulation/fragmentation 
model explained in Birnstiel et al. (2010a). The dust evolution 
is described by the advection-diffusion differential equation of 
the dust surface density S^, which for a single dust size can be 
written in cylindrical coordinates as 



dt 



]_d_ 
r dr 



(rS^Mr,d) ■ 



]_d_( y^j^d_ 
r dr\ * dr 



0, 



(1) 



where Dd is the dust diffusivity and Yg is the gas surface density. 
Because the timescales for gas viscous evolution are longer than 
the dust growth timescales, we considered that for the dust evolu- 
tion models, the gas surface density remains constant with time. 
This equation was solved for each size using the flux-conserving 
donor-cell scheme (see Birnstiel et al. 2010a, Appendix A). The 
radial velocity of the dust Mid has two contributions: the first one 
Mdrag comes from the drag with the gas that depends on the radial 
gas velocity Mr g and on the size of the particles; and the second 
one arises from radial drift Mdritt, which is proportional to the 
radial pressure gradient d,P, such that 



1 



Mr,d 



drP 



l+St^ St-i+StPgO 



(2) 



where the coupling constant St -the Stokes number- is defined 
as the ratio between the largest eddy turn-over time (I/O, with 
Q = ^jG Mi, r^^) and the stopping time of the particle within 
the gas. In the Epstein regime, where the ratio between the mean 
free path of the gas molecules /l^fp and the sizes of the particles 
a is /Imfp/fl > 4/9, St is defined at the disk midplane as 



St: 



aps n 



(3) 



where p, is the volume density of a dust grain, usually of about 
~ 1 g cm-^. 

In addition, the turbulent gas viscosity is considered as 
Y - Q'tmb c] ^ ' (Shakura & Sunyaev 1973), where c., is the 
isothermal sound speed, which is given by 



The interaction between the gas and the dust is fundamental for 
understanding how the particles evolve within the protoplanetary 
disk. When particles are well coupled to the gas, the dust relative 
velocities are mainly generated by Brownian motion and settling 
to the midplane. Considering these two sources for the velocities, 
particles stick by van der Waals forces and collisional growth 
is very efficient, producing mm-size grains in the outer regions 
(> 50 AU) on timescales of ~ 10^ years in typical T-Tauri disks 
(Birnstiel et al. 2010a). Unfortunately, when particles grow, the 



c? = 



keTir) 

IMlp 



(4) 



where k^ is the Boltzman constant, nip the proton mass and yu 
the mean molecular mass. Magnetorotational instability (MRI) 
is the most likely source of turbulence in disks (e.g Johansen & 
Klahr 2005). MRI essentially depends on the disk temperature 
and penetration of cosmic, X- and UV-rays to the midplane. The 
dust diffusivity in Eq. 1 can be defined in terms of St when the 
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Fig. 1. Vertically integrated dust density distribution (Eq. 6) after different times of evolution and including radial drift. Case of the 
BDdiskp-Ophl02pai-ameters (Mbo = 0.05 Mq, Led = 0.03 L^ and Tbd = 2880 /T) and 7?out = 30 AU, S = So r-°•^ 
Mdisk = 2 Mjup, v/ - 10 m i"' and cKturb = 10"^. The solid white line represents the particle size corresponding to St = 1 (Eq. 3) 
and reflects the shape of the gas density. The cyan line corresponds to the largest size that particles can reach given a fragmentation 
velocity v/. 



gas diffusivity is considered to be the turbulent gas viscosity v 
(Youdin & Lithwick 2007), hence 



Dh = 



1+St 



2' 



(5) 



The Stokes number St describes how well the particles are 
coupled to the gas. When St <K 1, the first term dominates in 
Eq 2, i.e., Wdiag' ^"d as a result these particles move along with 
the gas, meaning that they have sub-Kepleiian velocities. When 
the particle size increases, the second term, i.e., Mdrift starts to 
dominate and reaches the highest value when St = 1, and as a 
consequence these are the particles that react strongest to the 
sub-Keplerian velocity of the gas. The bodies with St » 1 move 
with their own velocity, i.e., Keplerian speed. 

In addition, dust particles grow, fragment, and crater de- 
pending on the relative velocities between them. For this the 
Smoluchowski coagulation equation (Smoluchowski 1916) was 
solved for the dust grain distribution n{a, r, z) (see Birnstiel et 
al. 2010a, Eq. 35 and Eq. 36), considering three different phys- 
ical processes: coagulation, fragmentation, and erosion. For the 
relative velocities, we assumed Brownian motion, settling, tur- 
bulent motion (Ormel & Cuzzi 2007), and drift velocities in the 
azimuthal and radial direction. The fragmentation velocities Vf 
were estimated based on laboratory experiments and theoretical 
work of collisions, which are of the order of few m s ' for sil- 
icates (Blum & Wurm 2008) and several m s"' for ices (e. g., 
Wada et al. 2009, 2011). To describe the dust grain distribution, 
we refer to the vertically integrated dust surface density distribu- 
tion per logarithmic bin, which is given by 



cr(r, a) 



Xoo 
CO 



(r, z,a) ■ m ■ adz. 



(6) 



where m is the mass of a single particle of size a. Therefore the 
total dust surface density can be written as 



^d(r) 



f 

Jo 



cr(r,a)dlna. 



(7) 



2.2. High radial drift and fragmentation in BD disks 

In a protoplanetary disk, large bodies feel no pressure and move 
with Keplerian speed, while the gas is slightly sub-Keplerian. As 
a result, particles that have grown to a size where they start to de- 
couple from the gas, i.e., St a; 1, feel a strong headwind, lose an- 
gular momentum, and move inward. Solving the radial Navier- 
Stokes equation, the azimuthal velocity of the gas is given by 
(Nakagawa et al. 1986) 



M^,g - Uk(l -1]), 

where m^ is the Keplerian speed and 

1 dP(r) 

2pg(r, z)rQ.^ dr 



(8) 



(9) 



with pgir, z) the gas density, such that the gas surface density is 

Xoo 
Pg{r, z)dz. For an ideal gas, the pressure P(r) 

is defined as 



P{r)^cipg{r,z), 
and therefore 

c] (d\np dlncl 

T] 1 

luk \dlnr dlnr 



(10) 



(11) 



where the term in parenthesis depends only on the exponents 
that characterize the power-law gas density and temperature ra- 
dial profiles in a flared disk. For the temperature, T{r) can be 
approximated by (Kenyon & Hartmann 1987) 



T(r) = T. 






1/4 



1/4 



-\fr 



(12) 



where r/Ri, > 10 . Cinc is the angle between the incident ra- 
diation and the local disk surface. For the last proportionality 
of Eq. 12 the weak dependence of ai,^^ with R^, is neglected. 
Considering Eq. 1 2 for the disk temperature and assuming a sim- 
ple power-law function for the gas density profile, the difference 



Pinilla P. et al.: Explaining millimeter- sized particles in brown dwarf disks 



between the orbital gas velocity and the Keplerian speed scales 
as (independently of the radial location r) 



^0,g ^k 



2uk 



dlnp c/lnCj 



dlnr dlnr 



T(r) 



-1/4 



(13) 



This difference determines how much angular momentum 
the particles lose, and as a consequence, how fast they drift to the 
central star. Taking Eq. 2 for the dust radial velocity, and consid- 
ering grains with St > 1, the drag term (first term of Eq. 2) can 
be neglected and Mi^ can be written in terms of the difference 
U(h o — iii_ as 



Mr d - 2 



St-i-HSt' 



(14) 



For comparison, taking the parameters of the BD know as 
-Oph 102 (Ricci et al. 2012) and comparing them with a 



young solar-type star, i.e., Mp_ophi02 = 0.05 



L'n-\ 



'p-OphI02 



= 0.03 X Lsun, this leads to 



(w^.g - Mi) 



p-OphI02 



2 I M0,e - Mi- 1 
V ^'^ /su 



and 



(15) 



The fact that the radial drift is higher for BD disks than for 
Sun-like disks also implies that destructive collisions due to ra- 
dial drift are more likely because it contributes to increase the 
relative velocities between the particles. Therefore rapid inward 
drift and the potential fragmentation due to these high drift ve- 
locities are major problems for particles in BD disks. 

In addition, the highest turbulent relative velocity of particles 
with a given Stokes number is given by (Ormel & Cuzzi 2007) 



/^ui 



3 fft, 



urb 



St + St 



-1 ■'■ 



(16) 



When the radial drift is set to zero, the largest size of particles 
«max is calculated when the turbulent relative velocities (Eq. 1 6) 
are equal to the fragmentation velocities vj. For particles with 
St < 1, Omax is approximately given by 



22, 



'f 



3;raturbPi c] 



(17) 



This implies that considering only turbulent velocities, a^a^ 
will depend on disk parameters such as o-turb, ^g, and T{r), which 
strongly vary between BD and T-Tauri disks. Whether turbu- 
lence or radial drift is the cause for destructive collisions for the 
dust within the disk, fragmentation would occur differently for 
BD disks, as we discuss in the following sections. 

As a general illustration of this problem in BD disks. Fig. 1 
shows the vertically integrated dust density distribution (Eq. 6) 
after different times of evolution, including radial drift and us- 
ing the dust evolution model described in Sect. 2.1. For this 
case, taking the parameters p-Oph 102 (Mbd - 0.05 Mq, 
Lbd - 0.03 Lq and Tbd - 2880 K), a truncated power 
law is taken for the gas density profile with 1. - Y.q r "^ and 
r e (1,30) AU such that the mass of the disk is Mdisk = 2 Mjup. 
The largest size that particles can reach is represented by the 
solid cyan line, which was computed considering the fragmenta- 
tion velocity as vj - 10 m 5"' and the turbulence parameter as 
ffturb - 10"^. The solid white line corresponds to St = 1, which 
reflects the shape of the gas surface density based on Eq. 3 and 
particles that feel the strongest radial drift. The dust particles 
initially grow in the disk, allowing to have mm size grains in the 



outer regions (r > 10 AU) after 0. 1 Myr of evolution (left panel 
of Fig. 1). However, when particles grow to sizes close to St ~ 1, 
turbulence and radial drift lead to fragmentation and inward mi- 
gration of particles, and in only 0.5 Myr of evolution (middle 
panel of Fig. 1), the outer region is empty of mm-grains. This 
scenario does not change significantly with time (right panel of 
Fig. 1), but the inner region (r < 10 AU) is even more depleted 
of mm grains, since they continue fragmenting and drifting to- 
ward the star, leading to a dust-poor disk after 1.5 Myr As a 
result, explaining how micron-size dust grows to pebbles and 
how those are retained in the outer regions of BD disks is very 
challenging. 

In addition, if perfect sticking and micron-sized compact 
particles are considered to be formed, the mean growth time 
of monomers with mass m at a fixed distance from the star is 
(Brauer et al. 2008; Okuzumi et al. 201 1) 



1 dm 
m dt 



(18) 



where dm/dt = pjacoiiAv, with pj is the dust density, cTcon 
the collision cross-section and Av the collision velocity. For a 
narrow-size distribution, perfect mixing and settling assumed 
in the Epstein regime, a given dust-to-gas ratio, the growth 
timescale is proportional to 






«Q-', 



(19) 



with Zrf as the dust surface density. The growth timescale in 
disks of BD and T-Tauri disks can be compared considering, for 
example, regions of similar temperature. Following Mulders & 
Dominik (2012), we re-scaled the distances with the stellar lu- 
minosity as well as the Keplerian frequency, to obtain 



grow y BD grow 



Led 
i7' 



(20) 



where t^^^ and tJ^^^ correspond to the growth timescales in 
BD and T-Tauri disks respectively. This implies that for micron- 
size particles and considering for simplicity only settling mo- 
tion, the mean collision time Tgjovv is generally shorter for BD 
than for T-Tauri disks. As an illustration, if we again take the 
BD p-Oph 102, the growth timescale due to settling is approxi- 
mately one order of magnitude shorter for dust particles around 
this BD than for particles in a Sun-like disk. As a consequence 
of the result that settling and radial drift occur faster in BD disks, 
the very early stages of dust growth, when particles just stick and 
grow through molecular interactions, is more efficient at the lo- 
cation of a given temperature in BD disks than in T-Tauri disks. 
Based on the low spectral indices (ai-smm < 3) measure with 
recent ALMA-Cycle 0-observations at 3.2 mm and 0.89 mm of 
the young BD p-Oph 102 (Ricci et al. 2012) and CARMA ob- 
servations of the BD disk 2M0444-I-2512 at 3 mm (Ricci et al. 
2013) and 0.850 mm (Bouy et al. 2008), we aim to explain how 
mm-grains can form and be retained in disks around BD. We 
divide this work in two parts: In the first part, we do not allow 
particles to drift to study whether we can create mm-grains con- 
sidering only grain growth and fragmentation. In the second part, 
we aim to explain how these grains are trapped in the outer re- 
gions of BD disks when radial drift is included. 
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Table 1. Model parameters 



Parameter 


Values 


MbdW^] 


0.05 


^bd[^g] 


0.03 


Ted [K] 


2880 


Pj [g/cm^] 


1.2 


Mdisk[Mji,p] 


2.0 


^i„[AU] 


1.0 


^ou,[AU] 


{15,30,60,1001 


Clurb 


{io-\io-*, 10-3) 


Vf [m/s] 


{10,30) 


A 


{0.4,0.61 


f 


1 



2.3. Set-up 

For the simulations of this work, we considered two different 
scenarios. In the first case, there is no dust radial evolution, i.e., 
the drag and drift terms in Eq. 2 are neglected. In the second 
case, we considered dust radial evolution with pressure bumps. 
For the gas density, we considered truncated power-law func- 
tions 2^ = So r"'', with /:> = (0.0, 0.5, 1.0) and Sq was com- 
puted such that the mass of the disk was always Mdisk = 2 Mjup. 
For the outer radius of the disk, four possibilities were consid- 
ered: 15, 30, 60 and 100 AU, since the exact typical disk spatial 
extent is unknown (Luhman et al. 2007; Ricci et al. 2012). 

The turbulence parameter aturb was taken to have values be- 
tween 10"^ - 10"^. This parameter directly influences the max- 
imum size of particles (Eq. 17), if turbulence is the cause of 
fragmentation. Mulders & Dominik (2012) inferred Cturb ~ 10""* 
from SED modeling, assuming a fixed grain size distribution and 
gas-to-dust ratio. The fragmentation velocities were assumed to 
be Vf = {10, 30) m s"', which are the values expected for ices 
(Wada et al. 2009). For all simulations, the dust-to-gas ratio was 
initially considered to be 1%, the initial size of the particles was 
taken to be 1 fim and the maximum size that particles can reach 
in the simulations was fixed to be 1 km, because for these sizes 
gravitational effects start to play a role, which are not included 
in the present dust evolution model. 

When we considered pressure inhomogeneities, we assumed 
the same prescription as in Pinilla et al. (2012a) to simulate long- 
lived pressure bumps, with the unperturbed density as a simple 
power law 



5:;(r) = So r 



l + A COS 



2;r- 



L(r) 



(21) 



where the wavelength of the perturbation L{r) is taken to be a 
factor / of the vertical disk scale-height h(r) i.e. L{r) - f h{r), 
with 



-1/8 



h(r) = Cj Q oc 



yfr. 



(22) 



Pinilla et al. (2012a) showed that under these assumptions, 
bumps with an amplitude of ~ 30% and a wavelength of one 
scale-height are necessary for T-Tauri disks, to reduce the radial 
drift and keep the particles in the outer regions of the disk after 
several Myr. However, in Sect 2.2, we showed that radial drift ef- 
fects are stronger for BD disks, independently of the distance of 
the particles from the star In addition, since the wavelength of 
the perturbation was considered proportional to the disk scale- 
height hi/), at a given distance r, L{r) would by definition be 



larger for BD than for T-Tauri disks (Eq. 22). This implies that 
for a given amplitude of the perturbed density, the pressure gra- 
dient would be lower in BD disks, making the trapping of par- 
ticles even more difficult. We used a fixed bump width equal to 
one scale height, i.e., / = 1 .0, since the scale of turbulent struc- 
tures from MRI are of the order of the scale height of the disk 
(Flock et al. 2011), and they may be the origin of pressure inho- 
mogeneities in the disk. The amplitudes considered are therefore 
higher than for T-Tauri disks and these are taken A={0.4,0.6). 
All stellar and disk parameters are summarized in Table 1 . 



3. Results 

In this section, we present the results of the numerical simula- 
tions setting the radial drift to zero, followed by the case of non- 
zero radial drift and pressure bumps in the disk. For each case, 
the observational perspectives are also presented. 

3.1. No radial drift 

3.1.1. Dust density distribution 

Figures 2 and 3 show the vertically integrated dust density distri- 
bution (Eq. 6) after 1 Myr of dust evolution when the radial drift 
is set to zero, for the case of Rout - 30 AU and two different 
values of the fragmentation velocity: 10 m s^^ and 30 m ,?"'. In 
each case, the three different values of the gas density slope p 
and the turbulence parameter o-turb are plotted. For this case, the 
dust is considered to be in a steady state, therefore Mdrag and Mdrift 
are neglected in Eq. 2. 

Effect of turbulence: when particles grow and reach sizes 
such that they are more weakly coupled to the gas, the main 
sources of relative velocities are turbulent and azimuthal veloci- 
ties, since radial drift velocities are set to zero. Because the high- 
est turbulent relative velocities (Eq. 1 6) depend on aturb and the 
disk temperature is assumed to be very low in BD disks (~ lOK 
at 10 AU), fragmentation is not caused by turbulent motions. 
Instead, it mainly happens because the azimuthal dust veloci- 
ties Ud^if, are as high as the limit beyond which particles fragment 
(fragmentation velocity, v/). The azimuthal drift velocity m^,^ is 
given by (Birnstiel et al. 2010a) 



Ud,il, — 



d,P I 1 

pfi 1 1 H- St- 



(23) 



Hence, as we noticed in Fig 2, the fragmentation limit and 
consequently a^^^^ are independent of the a^^th and for a given 
gas density slope, the largest size that particles can reach is the 
same for each atuib and only depends on the gas density profile. 
On the other hand, we notice that the concentration of mm- and 
cm-size grains is more evenly distributed in the whole disk for 
Ctui-b = {10^^, lO""*) than for aturb = 10"^, for which they are 
mainly in the inner part of the disk. 

The first remarkable difference of grain growth between T- 
Tauri and BD disks is that the turbulent mixing strength fftmb 
does not play as important a role for dust fragmentation in BD 
disks as it does for T-Tauri disks. Without radial drift, fragmen- 
tation in T-Tauri disks mainly occurs because of turbulent mo- 
tion and the effect of different aturb values is significant (see e.g 
Birnstiel et al. 2010b, for the effect of turbulence in the spectral 
index), while for BD disks the azimuthal relative velocity domi- 
nates, which is independent of fftuib- 
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Fig. 2. Vertically integrated dust density distribution after 1 Myr of evolution without radial drift. Case of Rom - 30 AU and 
Vf = 10 m s"' and different values of the gas density slope /:> = (top panels), p = 0.5 (middle panels) and p = 1.0 (bottom 

pane/ij; with the turbulence parameter atuib = 10^^ (left panels), atmb = 10^"^ (center panels) and axurb = 10^^ (right panels). 
The solid white line represents the particle size corresponding to St = 1 (Eq. 3) and reflects the shape of the gas density, while the 
cyan line represents the maximum size that particles can reach given a fragmentation velocity Vf. 



Effect of the gas density slope: For the same turbulence pa- 
rameter and different gas density slope p, the dust density dis- 
tribution changes (see Figs. 2 and 3). When p - I (steep sur- 
face density), large particles (a > 1 cm) are well coupled to 
the gas in the dense inner region and grow to even larger sizes 
(a ~ 50 cm) before they fragment. Therefore, with p - \ 
the mm-size particles are less concentrated in the region where 
the gas density is lower, i.e., in the outer region of the disk 
(r > 10 AU), than for p = {0.5, 1.0). 



Effect of the fragmentation velocity: when the velocity at 
which particles fragment increases to 30 m 5"' (Fig. 3), grains 
can grow to very large sizes such that St > 1 . Fragmentation does 
not occur and this implies that bodies that are no longer affected 
by gas (a > 1 m) grow and even in some cases reach the maxi- 
mum size of particles considered in these numerical simulations 
(~ 1 km). Small particles are not replenished and as a result the 
disk can be almost empty of small grains (< 1 cm). For exam- 
ple, when p - and cctuib — 10"^, most of the disk consists 
only of big grains (a > 10 cm). However, when p increases and 
the gas density is higher in the inner regions, larger particles are 
coupled to the gas and can still be affected by gas turbulence. 
For this reason, there is some mm-dust in the inner part for the 
same atui-b = 10"^ and/? - 1. 



As we discussed above, because of the low disk tempera- 
tures, destructive collisions due to turbulence are less likely in 
BD disks, and even more with these high fragmentation veloc- 
ities considered for ices (such as Vf - 30 m ,?"'). With ices, 
the material properties vary, leading to fragmentation velocities 
of around ~ 10-50 m s'^ (Wada et al. 2009, 201 1), but it is 
still a matter of debate if ices can have these high fragmentation 
velocities. Some scenarios produce disks with only large grains 
(Fig. 3), which may not allow a good agreement between these 
theoretical models and the recent submillimeter observations of 
disk around two BD. For this reason, we focus the following re- 
sults on the case of vy = 10 m ,?"', where the fragmentation 
allows a constant replenishment of small particles. 

Effect of the outer radius: When the outer radius of the disk 
decreases or increases, the same amount of dust is distributed 
in a smaller or larger region with the same Mjisk and Eg profile, 
and the possibility of having cm- mm-size particles in the outer 
regions changes, we discuss this effect below. 



3.1.2. Comparison with observations 

To compare model predictions with disk fluxes of low-mass 
disks, we considered the dust distribution from the simulations 
after they almost reached a quasi-static state ~ 1 Myr. We calcu- 
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Fig. 3. Vertically integrated dust density distribution after 1 Myr of evolution without radial drift. Case of Rom - 30 AU and 
Vf = 30 m s"' and different values of the gas density slope /:> = (top panels), p = 0.5 (middle panels) and p = 1.0 (bottom 

paneZij; with the turbulence parameter atuib - 10^^ (left panels), aturb = 10^"^ (center panels) and aturb = 10^^ (right panels). 
The solid white line represents the particle size corresponding to St = 1 (Eq. 3) and reflects the shape of the gas density. 



lated the opacities Ky for each grain size and at a given frequency 
V. We used for simplicity optical constants for magnesium-iron 
silicates (Jaeger et al. 1994; Dorschner et al. 1995) from the Jena 
database' and followed the Mie theory. Once the opacities were 
calculated, the optical depth Ty was computed as 



o-(r, a)Ky 



(24) 



where ; is the disk inclination, which we took to be zero. The 
flux of the disk at a given frequency Fy is therefore 



Fy^ 



2'^'^°^' r'fi.(r(.))[i -.-]... 



d^ 



(25) 



with d being the distance to the source, which was taken to be 
140 pc as in the young disks in TauiTis and Ophiuchus star- 
forming regions. By(T{r)) is the Planck function for a given tem- 
perature profile T{r). At mm wavelengths, the flux is propor- 
tional to the dust mass in the outer region of the disks. Hence, 
the flux Fy could be approximated as a power law Fy oc -/""<« ^ 
where the spectral index Cmm gives information about the size 
of the grains and is expected to be lower than 3 when the dust 
reaches mm sizes (Natta et al. 2007). Figure 4 shows the pre- 
dicted fluxes at 1 mm (Fimm) and the spectral index between 1 



and 3 mm (ffi-amm) for v/ = 10 m 5 
considered in Sect. 2.3. 



and all other parameters 



http://www.astro.uni-jena.de/Laboratory/Database/databases.html 



First of all, comparing Fimm and ai-smm for different Rom, 
we notice that the spectral slope increases when the disk is more 
extended. This is a natural result of distributing initially the same 
amount of dust in an extended region, decreasing the possibil- 
ity of having mm-size grains. The obtained fluxes are between 
~ 1-20 mJy for all cases. In addition, a single case (the only 
diamond-point of Fig. 4) where fragmentation does not happen 
is plotted, with v/ - 30 m s'\ Rout - 30 AU and p - 0.0. 
It is possible to see that most of the dust in the disk has grown 
to large sizes (a > 10 cm). When the particles have such a large 
size, they have very low opacities, which results in low millime- 
ter fluxes. The predicted fluxes are very low (close to 1 mJy) 
compared with the other fluxes (see top-left panel of Fig. 3 and 
the corresponding diamond-point in the second-left panel from 
the top to the bottom of Fig. 4). 

In general, when fragmentation still occurs, i.e., for 
Vf - lOm s"', the integrated flux weakly depends on the gas 
density slope p, therefore we focus in the subsequent results on 
the intermediate value, i. e., p - 0.5. In most of the cases, the 
spectral slope and the flux are slightly sensitive to the turbulence 
parameter Ctui-b, and ci-smm increases for low turbulence. This 
is because there are fewer mm-size grains that contribute to de- 
crease the spectral index when Ctuib is low (Fig. 2 and Fig. 3), as 
was discussed in Sect. 3.1.1. 
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Fig. 4. Predicted fluxes at 1 mm (/^imm) and the spectral index between 1 and 3 mm (ffi-smm) after 1 Myr of dust evolution, without 
radial drift, Vf - 10 m ,s ' and for all other parameters discussed in Sect. 2.3. Black dots with error bars are millimeter observations 
of the young BDp-Oph 102 (Ricci et al. 2012) and 2M0444+2512 (Ricci et al. 2013) 



Comparing the millimeter observations (dot points with er- 
ror bars in Fig 4) of two BD disks (2M0444+2512 and p 
Oph 102), it is important to notice that when the radial drift 
is set to zero, the predicted and the observed values of Fimm 
and Q;i_3nim match well for the brightest BD disk and in par- 
ticular for Rout = 15 AU. The error bars come from the 
the optical depth uncertainties for different grain composition 
(Beckwith et al. 1990). When the extension of the disk increases 
to Rout = 30 AU, the spectral slope has values that agree with 
observations (especially for higher values of o-turb), but the fluxes 
are slightly higher than those detected by mm-observations. This 
leads to the conclusion that in this case, a combination of frag- 
mentation with a minor drift is necessary to reduce the number of 
mm-grains and have a better agreement between theoretical pre- 
dictions and recent millimeter observations of 2M0444+2512. 
For the faintest disk p Oph 102, radial drift is indeed needed 
in any case. For /?out = 60 AU, only few cases allow to 
have ai-Bmrn ^ 3, nevertheless, these cases match the mil- 
limeter observations less well . Finally, for most of the cases of 
^out - 100 AU, the model predictions for the millimeter fluxes 
and spectral indices are inadequate to explain the observations. 

As a conclusion, theoretical models of dust evolution in 
which the radial drift is set to zero and considering BD disk con- 
ditions such as low-mass disks (M^isk - 2 Mj^p), low radial 
extension (/?out = 15 AU), presence of ices that allow fragmen- 
tation velocities of about v/ = 10 m .?"' and average turbu- 



lence strength Cturb = {10 "*, 10 ^), are the models with the best 
agreement to mm-observations. 



3.2. Radial drift and pressure bumps 

We focus our attention on the most favorable cases: 
Rout = {15,30,60) AU (particularly 15 and 30 AU) and 
ffturb = 10"^. The gas density slope was taken to be /7 = 0.5. 
The amplitudes for the long-hved pressure bumps were consid- 
ered A = {0.4, 0.6) (Eq. 21). Pinifla et al. (2012a) showed that 
these inhomogeneities are comparable with global simulations 
of zonal flows with an amplitude of ~ 25% by Uribe et al. (201 1). 
For the simulations with radial drift, M^jag and Mdnft for Eq. 2 
were taken into account. This implies that radial drift also con- 
tributes to the total relative velocities of dust particles. 



3.2.1 . Dust density distribution 

In a region where the pressure gradient is positive, gas moves 
with super-Keplerian velocity and particles with sizes corre- 
sponding to St ~ 1 would move outward. As a consequence, 
inside the pressure bumps the amount of dust naturally increases, 
allowing an increase in the frequency of sticking collisions. 
Figure 5 compares the dust density distribution after 1 Myr of 
evolution for the case of /?out = 30 AU, v/ = 10 m s"'. 
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Fig. 5. Vertically integrated dust density distribution after 5 Myr of evolution including radial drift and a bumpy gas surface density 
(Eq21)witii A- OA (left panel) and A = 0.6 (right panel). Case of Ro^^ = 30 AU, vy = 10 ms"',;? = 0.5 and aturb = lO"-'. 
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cyan line corresponds to the fragmentation Umit. 
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Fig. 6. Dust-to-gas ratio at different times of dust evolution in- 
cluding radial drift and a bumpy gas surface density (Eq 2 1 ) with 
A = 0.4. Case of atui-b = 10"^^'/ = 10 m s"', p = 0.5, 
^out — 30 AU (top panel) and Ro^t = 60 P^J (bottom panel). 



p = 0.5, atuib - 10"^ and two different values of the amplitude 
of the perturbation A = 0.4 and A = 0.6. In both cases, with 
the pressure bumps considered, it is possible to reduce the strong 
radial drift and retain mm-size particles even after million-year 
timescales. The efficiency of the pressure bumps is evident, with 
many more mm grains than in the test case of Fig. 1 (empty 
disk after 0.5 Myr of dust evolution. Fig. 1). In addition, the ef- 
fectiveness of the amplitude strength, i.e., the pressure gradient 
strength, for the trapping of particles is almost the same for both 
amplitudes after 1 My of evolution. 

Inside the bumps, the radial drift is reduced and fragmen- 
tation is also less efficient than in the case of Fig. 1. With a 



fragmentation velocity of vy = 10 m s , particles grow to 
sizes corresponding to St ~ 1 for the cases in Fig. 5, and be- 
cause of the positive pressure gradient, grains are trapped inside 
the bumps. However, since fragmentation due to turbulent and 
azimuthal relative velocities is still happening, when particles 
collide and become smaller such that St < 1, they are more diffi- 
cult to trap in the pressure bump (Pinilla et al. 2012a) and would 
finally drift to the star 

3.2.2. Dust-to-gas ratio 

Figure 6 shows the radial dependence of the dust-to-gas ratio 
at different times of evolution for the cases of atuib = 10"^, 
v/ = 10ms"',andp = 0.5, in a bumpy gas surface density 
with A = 0.4 and Roni = {30,60) AU. For this case, the 
dust-to-gas ratio remains almost constant with time until 1 Myr 
of evolution. It increases inside the pressure bumps and varies 
between ~ 10"^ - 10"' , but there is a remarkable reduction after 
5 Myr of evolution, where it can reach values of 10"^ in the inner 
part and 10 "* - 10"^ in the outer parts of the disk. 

When the outer radius increases, the dust-to-gas ratio does 
not remain constant with time in the outer regions of the disk, 
and it is possible to distinguish that it decreases with time, reach- 
ing even values of ~ 10"* in those regions (bottom panel of 
Fig. 6). 

It is important to note that the gas density profile remains 
constant for our dust evolution models, since viscous evolution 
timescales are longer than dust growth timescales. In addition, 
we did not consider any mechanism that might selectively dis- 
perse the gas, such as photoevaporation (see e.g. Gorti et al. 
2009; Owen etal. 2011). 

3.2.3. Comparison witli millimeter observations 

Focusing the attention on the cases where fragmentation is ef- 
fective (i.e., when v/ = 10 m ,?"' and ofturb - lO""^), Fig. 7 
shows the time evolution of Fimm and ffi-amm including radial 
drift and a bumpy gas surface density with 7?out - {15,30)AU, 
p = 0.5, and A = {0.4,0.6). For A = 0.4, we notice a good 
match between the theoretical predictions and the observations, 
even after 1 Myr of evolution, there are still enough mm-grains 
in the outer regions that contribute to have high values of the flux 
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at 1 mm. The fluxes become lower at 5 Myr, because a signifi- 
cant amount of grains continues to fragment and drift toward the 
star. However, Fimm and ffi-smm at intermediate time steps like 
1-2 Myr coincide well with the observed values, and are in the 
range of the estimated ages of the observed BD disks. The ef- 
fect of increasing the pressure gradient, i.e., higher A, is notable, 
since the fluxes increase because there are more trapped grains, 
in particular at long timescales, for instance 5 Myr. 

To achieve a better agreement of the predicted fluxes Fimm 
and the spectral index ai-smm with the current mm-observations 
of 2M0444-I-2512 and p Oph 102, it would be necessary to 
increase the mass of the disk. Figure 8 shows the compar- 
ison of Fimm and the ai-imm for two different disk masses 
(A^disk = {2,5) Mjup) and the specific case of /?out - 15 AU, 
p = 0.5, O'tui-b - 10^1^/ - 10 mi"' and A = 0.4. Because 
there is more dust for Mdisk - 5 Mjup than for Mdisk = 2 Mjup, 
there are more mm-grains that contribute to the integrated flux, 
and therefore the spectral slope can be lower. 

4. Summary and discussion 

One of the most important problems in the core accretion the- 
ory for planet formation is "the meter-size barrier". This phe- 
nomenon results from the combination of high radial drift and 
fragmentation. A meter-sized object at 1 AU drifts toward the 
star on timescales much shorter than the growth timescales. In 
the colder regions of a disk, the same problem occurs when 
translated to mm-sized grains. Both radial drift and fragmenta- 
tion appear to be different for particles in BD disks. 

We have studied how dust evolves and explored diff'erent sce- 
narios where mm-grains survive in BD disks. For the results pre- 



sented here, we considered typical BD parameters. For instance, 
the disk outer regions reach temperatures of ~ 10 K. Assuming 
that MRI is the source of turbulence, which depends basically 
on the disk temperature and considering o-tuib-prescription, low 
and moderate values of fftuit were taken (10^^ - lO""*). In ad- 
dition, the presence of ices in the disk was assumed and frag- 
mentation velocities were taken to be v/ = {10,30) m s^\ 
based on numerical simulations of collisions with ices (see e.g. 
Wada et al. 2009). Furthermore, we considered low-mass disks 
(2 Mjup) and typical values for masses and luminosities of BD 
(see Table 1). For the radial extension of BD disks, we consid- 
ered /?out = {15,30,60, 100) AU. Luhman et al. (2007) pre- 
sented observations combining Spitzer spectroscopy and high- 
resolution imaging from HST of a circumstellar disk that is in- 
clined close to edge-on around a young brown dwarf in Taurus, 
estimating a disk radius of ~ 20-40 AU. ALMA-Cycle obser- 
vations of p-Oph 102 (Ricci et al. 2012) do not resolve the disk, 
but they estimate an outer radius for the dust of ~ 15-40 AU. 
From 3-mm CARMA observations of 2M0444-I-2512 (Ricci et 
al. 2013), the outer radius for the dust is estimated to be ~ 20- 
75AUforp = 1 and ~ 15-30 AU for p = 0.5. 



4. 1 . Fragmentation and drift barriers 

We have shown how several dust growth aspects in BD disks 
differ from T-Tauii and Herbig Ae/Be disks. In the first stages 
of dust evolution, when particles move along with the gas and 
coagulate, relative velocities are mainly due to Brownian motion 
and settling to the midplane. It was discussed in Sect. 2.2 that 
growth timescales due to settling are shorter at the location of a 
given temperature in BD than in T-Tauri disks. 



10 
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When particles reach sizes such that they start to decouple 
from the gas (i.e. St ~ 1), radial drift and turbulent motion be- 
come the main sources for relative velocities. On one hand, we 
demonstrated in Sect. 2.2 that radial drift can be even twice as 
fast for particles in BD disks than in T-Tauri disks, and as a con- 
sequence they would be depleted on shorter timescales. On the 
other hand, considering moderate Ctuib and low disk tempera- 
tures, turbulent relative velocities are quite low. If turbulence is 
the main cause for fragmentation, destructive collisions are less 
likely in BD disks than in T-Tauri disks. When radial drift and 
fragmentation are considered, BD disks would be dust-poor af- 
ter only few Myr because of the fast inward migration of dust 
grains (Fig. 1). 

4.2. No radial drift 

This drastic case, where radial drift is not included, is important 
for studying whether mm-grain can be indeed formed consid- 
ering dust coagulation/fragmentation models. When radial drift 
was set to zero, we found that two different scenarios can hap- 
pen under BD disk conditions. In one scenario fragmentation 
occurs because of the relative azimuthal velocities, considering 
Vf = 10 m i"' (Fig. 2). In this case, it is possible to form 
mm-size particles in BD disk and achieve a good match for the 
predicted mm-fluxes and spectral indices, in particular for the 
brightest BD disk 2M0444H-2512 (Fig. 4) for R^ut = 15 AU. 
However, for larger disk sizes, i.e., 7?out = {30, 60) AU, a slight 
radial drift is necessary to decrease the predicted mm-fluxes. For 
the other BD diskp Oph 102, radial drift is always needed to re- 
duce the mm-fluxes by at least a factor of 3. 

The second scenario is when fragmentation is less likely, i.e., 
Vf = 30 m s^^ (see e.g. Fig. 3), when most of the disk consists 
of big grains (a > 10 cm), leading to low mm-fluxes that un- 
derestimate millimeter observations. Therefore, values such as 
Vf = 30 m ,?"' do not yield a good agreement between mod- 
els and observations, and average values of Vf - 10 m 5"' are 
needed. 



4.3. Radial drift and pressure bumps 

To reduce the rapid inward drift that dust particles experience in 
protoplanetary disks and explain the presence of mm-grains in 
the outer regions of disks, pressure bumps have been suggested 
as a possible solution (see e.g. Pinilla et al. 2012a). In this work, 
when radial drift was taken into account, we considered strong 
gas density inhomogeneities (40% and 60% of amplitude), that 
led to regions in the disk with a positive pressure gradient. These 
pressure bumps could be the result of MRI effects (Johansen et 
al. 2009; Uribe et al. 201 1). However, there are no known mech- 
anisms that could produce this type of pressure bumps with these 
high amplitudes globally present in the disk. Nevertheless, if the 
scale-height of a disk is higher, as in the case of BD disks, the 
scale of turbulent structures increases, leading to higher pressure 
scales (e.g. Flock et al. 2011). 

In addition to global pressure inhomogeneities from MRI in 
the disk, a single and high-pressure bump can exist in disks. For 
instance, Rossby wave instability (RWI) generated by the pres- 
ence of a dead zone (Regaly et al. 2012) would create a high 
single-pressure bump. The possibility of a dead zone in a BD 
disk is still debated. On the other hand, the presence of a mas- 
sive planet could be also the reason for trapping of particles in a 
single huge pressure bump (Pinilla et al. 2012b). However, from 
low-mass disks, it is unlikely to have planets massive enough to 
open a gap. The criterion to open a gap in a disk depends on the 
local scale height h(r), meaning that a clear gap would be less 
likely for BD than T-tauri disks for the same disk and planet pa- 
rameters (Crida et al. 2006). Hence, these two possibilities were 
ruled out for the conditions that we considered for BD disks, and 
we only focused on global-synthetic pressure bumps that allow 
a broad radial distribution of millimeter dust grains in the whole 
disk. 

For simulations with pressure bumps, we considered a 
sinusoidal perturbation for the gas surface density with two 
different amplitudes A = {0.4,0.6) (Eq 21) and a wavelength 
equal to one disk scale-height. The parameters considered for 
these simulations were taken from the most optimistic no-drift 
cases. The pressure bumps cause an accumulation of dust in 
the location of pressure maxima, where radial drift is reduced 
and fragmentation due to radial drift decreases. Inside those 
regions, fragmentation would happen only because turbulent 
or azimuthal motions and a good match between theoretical 
perspectives and spectral slopes found with millimeter ob- 
servations for the two BD disks is possible, especially for 
^out - {15,30) AU. Increasing the mass of the disk from 2 
to 5 Mjup, theoretical predictions and mm-observations match 
even better. The disk mass predictions for p-Oph 102 (Ricci et 
al. 2012) are 0.2-0 .6 Mjup and for 2M0444+2512 (Ricci et al. 
2013) are 2.0-5.0 Mjup, assuming a gas-to-dust mass ratio of 100. 



4.4. Comparison between T-Tauri and BD disks 

Although T-Tauri and BD disks have very different properties, 
a first-order comparison between the models which agree well 
with millimeter observations in each case can be made. If we 
compare the best model described in Pinilla et al. (2012a) for 
T-Tauri disks (pressure bumps with an amplitude A of 30% and 
one scale-height of the disk as the width / of the bumps) with 
the best model for BD disks (A = 0.4, / = 1.0, Rom = 15 AU and 
A^disk = 2 Mjup), considering the same disk and dust properties 
(fragmentation velocity (vy = 10m s"'), viscosity (ofturb - 10"^), 
initial dust-to-gas ratio(l/100), etc), we find that 85% of the dust 
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mass is in grains larger than 1 mm in the BD case, versus 70% 
for T-Tauri disks. As a consequence, we expect our models to 
produce slightly flatter SEDs for BD disks than for TTS disks. 
We stress that this conclusion, which is inferred from our 1-D 
models that are only relevant in the mid-plane, should be taken 
with care. A proper vertical structure calculation is needed to 
assess this question, and see if the effect is strong enough to ac- 
count for 15-20% variation reported by e.g. Mulders & Dominik 
(2012) on Chameleon I targets. 

4.5. Additional improvements 

The aim of this paper was to determine the best conditions for 
grain growth in the mid-plane of BD disks. Our 1-D model is 
well-suited to this goal, but requires to azimuthally average the 
surface density profiles, which implies that we do not take into 
account possible local disk features that may locally affect dust 
growth. 2-D and 3-D calculations would allow the vertical disk 
structure and vertical grain distribution to be computed and pro- 
vide us predictions at shorter wavelengths (e.g. shape of the SED 
and of the silicate feature) that can be directly compared with 
large surveys. 

Our model can also be improved in the prescription of the 
dust-evolution process. Recent laboratory experiments using sil- 
icates have revealed that particles probably also bounce (Guttler 
et al. 2010), which is not included in our model so far. It is an 
open question however whether this happens in the case of colli- 
sion between ices (Wada et al. 201 1), and whether it has a strong 
impact on dust growth. In fact, the bouncing barrier (Zsom et al. 
2010) can be overcome if Maxwellian velocity distributions of 
relative velocities among dust grains are considered (Windmark 
et al. 2012a). Windmark et al. (2012b) also showed that after 
including bouncing effects and considering mass-transfer colli- 
sions, km-sized objects can still form at few AU from the star, 
but only if radial drift is neglected and cm-sized seeds are in- 
serted initially. 

Finally, the main caveat of our study is that we assumed that 
the gas density profile remains static for million-year timescales, 
neglecting the impact of different physical processes, for exam- 
ple the potential time evolution of the pressure bumps, like the 
case of zonal flows (see e.g. Uribe et al. 2011). Simultaneous 
modeling of gas and dust evolution is a step forward of our mod- 
els. 



5. Conclusion 

The first steps of planet formation differ between disks around 
T-Tauri stars and disks around BD. In BD disks, settling is more 
efficient and they are expected to be flatter. In addition, when 
dust grows and particles start to decouple from the gas, radial 
drift is a major problem and fragmentation due to turbulence is 
less likely than in T-Tauri disks. Our models show that, when 
BD disks are small (e.g.. Rout S 15AU), radial drift needs to be 
completely suppressed to account the mm-observations of one of 
the two BD disks observed so far. When disks are larger, a small 
degree of radial drift is necessary for both cases. Including radial 
drift, the most favorable conditions after million-year timescales 
of dust evolution are with strong pressure inhomogeneities of 
an amplitude of ~ 40%-60%. Global dust evolution models that 
include the vertical structure are the next step to reproduce all 
observational signatures of BD disks, for example the flatter ge- 
ometry of BD compared with T-Tauri disks, or the large average 
grain sizes in BD disk atmospheres. 
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